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ADAPTIVE FINITE ELEMENT APPROXIMATION OF 
STEADY FLOWS OF INCOMPRESSIBLE FLUIDS WITH 
IMPLICIT POWER-LAW-LIKE RHEOLOGY 

CHRISTIAN KREUZER AND ENDRE SULI 


Abstract. We develop the a posteriori error analysis of finite element approx¬ 
imations of implicit power-law-like models for viscous incompressible fluids. 
The Cauchy stress and the symmetric part of the velocity gradient in the class 
of models under consideration are related by a, possibly multi-valued, maximal 
monotone r-graph, with < r < oo. We establish upper and lower bounds 
on the finite element residual, as well as the local stability of the error bound. 
We then consider an adaptive finite element approximation of the problem, 
and, under suitable assumptions, we show the weak convergence of the adap¬ 
tive algorithm to a weak solution of the boundary-value problem. The argu¬ 
ment is based on a variety of weak compactness techniques, including Chacon’s 
biting lemma and a finite element counterpart of the Acerbi—Fusco Lipschitz 
truncation of Sobolev functions, introduced by L. Diening, C. Kreuzer and E. 
Siili [Finite element approximation of steady flows of incompressible fluids with 
implicit power-law-like rheology. SIAM J. Numer. Anal., 51(2), 984-1015]. 


1. Introduction 

Typical physical models of fluid flow rely on the assumption that the Cauchy 
stress is an explicit function of the symmetric part of the velocity gradient of the 
fluid. This constitutive hypothesis then leads to the Navier-Stokes system and 
its nonlinear generalizations, such as various electro-rheological flow models. It is 
known however that the framework of classical continuum mechanics, built upon 
the notions of current and reference configuration and an explicit constitutive equa¬ 
tion for the Cauchy stress, is too narrow for the accurate description of inelastic 
behavior of solid-like materials or viscoelastic properties of materials. Our starting 
point in this paper is therefore a generalization of the classical framework of con¬ 
tinuum mechanics, referred to as implicit constitutive theory, which was proposed 
recently in a series of papers by Rajagopal and collaborators; see, for example, 
|Raj03[ [Raj06[ [RS08| . The underlying principle of implicit constitutive theory in 
the context of viscous flows is the following: instead of demanding that the Cauchy 
stress is an explicit function of the symmetric part of the velocity gradient, one 
may allow an implicit and not necessarily continuous relationship between these 
quantities. This then leads to a general theory, which admits fluid flow models with 
implicit and possibly discontinuous power-law-like rheology; see, |Main7l IMainS] , 
Very recently a rigorous mathematical existence theory was develo ped for these 
models by Bulicek, Gwiazda, Malek, and Swierczewska-Gwiazda in |BGMSG09 , 
for r > for the range < r ^ the Acerbi-Fusco Lipschitz truncation 
[AF88| was used in order to prove the existence of a weak solution. In |DKS13| . 
using a variety of weak compactness techniques, we showed that a subsequence of 
the sequence of finite element solutions converges weakly to a weak solution of the 
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problem as the finite element discretization parameter h tends to 0. A key new tech¬ 
nical tool in the analysis presented in |DKS13| was a finite element counterpart of 
the Acerbi-Fusco Lipschitz truncation of Sobolev functions. However, in the case of 
velocity approximations that are not exactly divergence-free the convergence theory 


< r < CO. 


developed there was restricted to the range 

The focus of the present paper is on the adaptive finite element approximation 
of implicitly constituted power-law-like models for viscous incompressible fluids. As 
in |DKS13| . the implicit constitutive relation between the stress and the symmetric 
part of the velocity gradient is approximated by an explicit (smooth) constitutive 
law. The resulting steady non-Newtonian flow problem is then discretized by a 
mixed finite element method. Guided by an a posteriori error analysis, we propose 
a numerical method with competing adaptive strategies for the mesh refinement 
and the approximation of the implicit constitutive law, and we present a rigorous 
convergence proof generalizing the ideas in |MSV07| and |Sie09| . More precisely, we 
show that a subsequence of the adaptively generated sequence of discrete approxi¬ 
mations converges, in the weak topology of the ambient function space, to a weak 
solution of the model when < r < oo. In contrast with |DKS13| . stimulated by 
ideas from |BGMSG12 we shall be able to avoid resorting to the theory of Young 
measures. We emphasize that even in the case when the weak solution is unique we 
have only weak convergence of a subsequence; in this case, however, such a subse¬ 
quence can be identified with the aid of the a posteriori bounds derived herein; cf. 
Remark [31] 

The paper is structured as follows. In Section [^ we shall formulate the prob¬ 
lem under consideration and will introduce some known mathematical results. In 
Section we define the finite element approximation of the problem and present 
related technical properties and tools, such as the discrete Lipschitz truncation from 
|DKS13| . Section is concerned with the a posteriori error analysis for both the 
error in the approximation of the graph and the finite element approximation. The 
adaptive algorithm together with our main result are stated in Section for the 
sake of clarity of the presentation certain technical parts of the proof are deferred to 
Section We conclude the paper by discussing concrete graph approximations for 
certain problems of practical relevance. While the emphasis in this paper is on the 
mathematical analysis of adaptive finite elements algorithms for implicitly consti¬ 
tuted fluid flow models, the ideas developed herein may be of more general interest 
in the convergence analysis of adaptive finite element methods for other nonlinear 
problems in continuum mechanics with possibly nonunique weak solutions. 


2. Implicitly constituted power-law-like fluids 

In this section we introduce the variational model of the steady flow for an 
incompressible fluid with an implicit constitutive law given by a maximal m onotone 
x-dependent r-graph. We then recall the existence result from [BGMSG09 together 
with some known results and mathematical tools from the literature. Finally we 
conclude the section by discussing the approximation of x-dependent r-graphs by 
regular single-valued tensor fields. 


2.1. The variational formulation. Before we are able to state the weak formula¬ 
tion of the problem we need to introduce basic notations and recall some well-known 
properties of Lebesgue and Sobolev function spaces. 

For a measurable subset ui c we denote the classical spaces of Lebesgue 
and vector-valued Sobolev functions by and s e 

[l,oo], respectively. We denote the closure of := in by 

and let R)}:= {v e divr; = 0}. Moreover, we denote 

the space of functions in L^{uj) with zero integral mean by Lq{uj). For s, s' e (1, oo) 
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with i -I- ^ = 1 we have that L® (17) and Lg (17) are the dual spaces of L®(17) and 
Lg(17), respectively. The dual of Wg^’^(17)‘^ is denoted by For w = 17 

we omit the domain in our notation for norms; e.g., we write ll-llj, instead of ||•||g q. 

Let 17 c d G N, be a bounded open Lipschitz domain with polyhedral bound¬ 
ary. For r G (1, oo), we define r' G (1, oo), with 7 + ^ = 1, and set 


r 



dr 

d—r 


if r ^ 
^ d+2 

otherwise. 


( 2 . 1 ) 


With such r, r' and r, we shall consider the following boundary-value problem. 
Problem. For f e U (17)^^ find {u,p, S) G VFg’’'(17)'^ x Lq(17) x U (17)'^^'^ such that 


div(M ®u + pi — S) = f 
div u = 0 

(DM(a;), S{x)) g A{x) 


in I?'(17), 
in X>'(17), 

for almost every a; G 17. 


( 2 . 2 ) 


Here, Dit := ^(Vtt + (Vit)"^) signifies the symmetric part of the gradient of u. As 
is indicated in our choice of the solution space for the velocity u in the statement 
of the above boundary-value problem, we shall suppose a homogenous Dirichlet 
boundary condition for u. 

The implicit constitutive law, which relates the shear rate to the shear stress, 
is given by an inhomogeneous maximal monotone r graph A : x A{x). In 
particular, we assume that for almost every a; G 17 the following properties hold: 

(Al) (0,0 )gA(x); 

(A2) For all (<5i, cTi), (^ 2 , cr2) G A(x), 

{(Ti — (T2) : (di — ^2) > 0 (- 4 (a;) is a monotone graph), 

and if ¥= S2 and er^ ¥= <T2, then the inequality is strict; 

(A3) If (d, cr) G X and 

{& — cr) : (d — d) ^ 0 for all (d, a) G A{x), 
then (S,cr) e A{x) (i.e., A{x) is a maximal monotone graph); 

(A4) There exists a nonnegative function m G L^(17) and a constant c > 0, such 
that for all {6, cr) g A{x) we have 

cr : 6 ^ —m{x) -f c(|d|’' -f |cr|’' ) (i.e., A{x) is an r-graph); 

(A5) The set-valued mapping A : 17 —> x is measurable, i.e., for any 

closed sets Ci,C 2 c we have that 

{t G 17 : A{x) n (Ci x C 2 ) A 0} 
is a Lebesgue measurable subset of 17. 


Remark 1. Two remarks on the definition of x-dependent maximal monotone graph 
A are now in order: 

• Let T c X be a closed set; then, for a.e. x e LI we have that the set 

A{x) nT is closed. To see this, we assume w.l.o.g. that A{x) nV ¥= 0 and let 
{(dfc,crfc)}fegN A{x) n V, such that Sk ^ S e and cr^ cr e as 

n —> cx). Let (S,&) G A(x) be arbitrary. We then have that 

0 sS (d- - (Tfe) : (d - 5k) {a- - cr) : {5 - 5) 

as k ^ CO. This proves that {5,cr) G A{x) thanks to condition and the 

closedness ofD. Taking T) := {cr} x we then deduce that the set 

{rTeRty^-.{5,a)eA{x)} 
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is closed. This is condition (A5)(i) of |BGMSG09| . 

• According to |AF091 Theorem 8.1.4] property \(A5)\ is equivalent to the fact that 
the graph of the set-valued map A{x) belongs to the product a-algebra £,(fl) (g) 
® Here £(11) denotes the Lebesgue measurable subsets of 

and 18 the Borel subsets With the same argument it follows 

that (A5) is equivalent to the fact that, for any closed C c the sets 

{(x, (t) e n X : there exists d e C, such that {S, er) e A{x)], 

{(cc, e H X : there exists cr e C, such that {5, cr) e A{x)'^ 

are measurable relative to £(H) (g) 18(Msy^). These equivalent conditions imply 
that there exist measurable functions (so-called selections) S*, D* : H x —> 
such that (<5, S^{x, d)),{D^ {x, cr), cr) e A{x) for a.e. a; e H and all S,cr e 
^symi compare also with |BGMSG12[ Remark 1.1]. 

Buh'cek, Gw iazda, Malek , and Swierczewska-Gwiazda proved the following exis¬ 
tence resnlt in |BGMSG09| . 

Proposition 2. For r > there exists some (not necessarily unique) weak solu¬ 
tion of (2.21. 


2.2. Analytical framework. We shall shortly recall some resnlts that are crucial 
for the existence theory of (2.21. 

Inf-sup condition. The inf-sup condition has a central role in the analysis of the 
Stokes problem. It states that, for s, s' e (l,oo) with --1-^ = 1, there exists an 
Os > 0 such that 

q div V dx 


sup 

(viy 


> as llgjlj,, for all qe Lq {Ft). 


Ill,s 


This follows from the existence of the Bogovskii operator 18 : To(^) 
with 


(2.3) 


<’^( 11 ), 


divl8/i = h and Os ||18/i||j^ ^ ^ \\h\\,. 

for all s e (l,oo); compare e.g. with [DRSlOl lBog79| . It follows from |BF91I §11, 
Proposition 1.2] that condition (2.3) is equivalent to the isomorphism 

L«'(fl) ^ |n' G (v', w) = 0 for all w G ■ (2-4) 


Korn’s inequality. According to (2.2) the maximal monotone graph defined in 


(Al) (A5) provides control only on the symmetric part of the velocity gradient. 


Korn’s inequality states that this is already enough to control the norm of a Sobolev 
function; i. e., for s G (1, oo), there exists a > 0 such that 


7s ||f Ill s < l|D'*^lls ^ (2.5) 

compare e. g. with |DRS10] . 

We conclude this subsection with Gha con’s biting lemma and a corollary of it 
that is relevant for our purposes |GMS07l Lemma 7.3]. 


Lemma 3 (Ghacon’s biting lemma). Let Ft be a bounded domain in and let 
{fn}neN be « bounded sequence in L^{Fl). Then, there exists a nonincreasing se¬ 
quence of measurable subsets Ej cz Ft with \Ej \ —> 0 as j —> oo, such that is 

precompact in the weak topology of L^{Fl\Ej), for each j G N. 

In other words, there exists a v e L^{Fl), such that for a subsequence (not rela¬ 
belled) of Vn ^ V weakly in L^{Fl\Ej) as n —» oo for all j G N. We denote 

this by writing 


Vn V in L^{Fl) 
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and call v the biting limit of {vn}neN- 

Lemma 4 f |GMS07[ Lemma 7.3]). Let {fn}neN L^{Ll) be a sequence of nonneg¬ 
ative functions such that Vn v for some v 6 L^{n). Then, 


lim [ {vn-v)dx = 0 implies that Vn ^ v weakly in L\n) as n ^ co. 

2.3. Approximation of maximal monotone r-graphs. In general an x-depen- 
dent maximal monotone r-graph A satisfying |(A1) (A5) cannot be represented in 
an explicit fashion. However, it can be approximated by a regular single-valued 
monotone tensor field based on a regularized measurable selection S* with the fol¬ 
lowing properties; compare with |BGMSG09[ IBGMSG12 and Remark 

Lemma 5 (' |BGMSG12l Lemma 2.2]). Let S* : x x ^ be 

a measura ble selection of the x-dependent maximal monotone r graph A with the 
Then, for S,cr e 
for almost all x e Ll: 


properties (A1)-(A5) 
equivalent 


the following two statements are 


• (cr - S*{x,D)) 

• {S, (t) 6 A{x). 


(S — D) ^ 0 for all D e 


^dxd 

^sym> 


However, our convergence analysis relies on stronger structural properties of the 
graph approximation. We shall demonstrate in Section]^ how such graph approx¬ 
imations can be constructed for some typical problems of practical interest within 
the class of problems under consideration. 

Assumption 6. For n e N, there exists a mapping S” : x ^ such 

that 

• : LI ^ Klsym is measurable for all S e Kgym; 

• S^{x, •) : Mgym ^ I^sym is continuous for almost every x e LI; 

• is strictly monotone; i.e., for almost every x e 12 we have 

(S”(x,di) - 5 "(x,< 52)) : (<5i - ^ 2 ) > 0 for all 61 ^ 62 ^ 

• There exist constants ci,C 2 > 0 and nonnegative functions rh e L^(L1), 
k e L^ {LI) such that, uniformly in n e N, we have 

|iS"(x, 5)1 ^ Cl 1^1^ ^+ k{x) and S^”(x, 5) : 5 > C 2 — to(x) 
for all 8 e and almost every x e 12. 

We aim to approximate the solution of \2.2\ by solving the following explicit 
nonlinear problem: For f e L^ (12)'^ find {u,p, S) e Wg^’’^(12)‘^ x L]j(12) x L'" (12)'^^'^ 
such that 

div{u(^u + pi — S) = f in I?'(12)‘^, 

divit = 0 in 22'(12), (2.6) 

S{x) = S'"(x,Dm(x)) for almost every x e 12. 


3. Finite Element Approximation 


This section is concerned with approximating problem (2.21 by finite elements. 
To this end, we first approximate (2.2) by an explicitly constituted problem. We 
introduce a general finite element framework for inf-sup stable Stokes elements. 
This, together with some representative examples of velocity-pressure pairs of finite 
element spaces, is the subject of section 13.2] The finite element approximation 


of (2.2) is stated in section 3.3 
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3.1. Domain partition and refinement framework. In this section we provide 
the framework for adaptive grid refinement. For the sake of simplicity, we restrict 
our presentation to conforming simplicial meshes and refinement by bisection. To be 
more precise, let Qq be a regular conforming partition of into closed simplexes, the 
so-called macro mesh. Each simplex in the partition is referred to as an element. We 
assume that there exists a refinement routine REFINE with the following properties. 

• The refinement routine has two input arguments: a regular conforming partition 
Q and a subset M cz Q of marked elements. The output is a refined regular 
conforming triangulation of fl, where all elements in Ai have been bisected at 
least once. The input grid can be Qo or the output of a previous application of 
REFINE. 

• Shape-regularity: We call Q' a refinement of Q (briefly Q' ^ Q), when it can be 
produced from ^ by a finite number of applications of REFINE. The set 

G ■.= {Q: Q is a refinement of ^ 0 } 

is shape-regular, i. e., for any element E e Q, Q e G, the ratio of its diameter to 
the diameter of the largest inscribed ball is bounded, uniformly with respect to 
all partitions Q, Q e G. 

For the existence of such a procedure, we refer to |Ban MM os94| or the monograph 
|SS05| and the references therein. 

The neighbourhood of an element E e Q, with ^ e G, is denoted by 
Af^{E) := {E' eg-. E' nE ^ 0}. 

Let a; c n and define (uj) := [J {E e Q \ E n uj ¥= 0}. For subsets M. ez Q, let 

n{M):=[j{E\EeM} and (M) := {n{M)). 

Thanks to the shape regularity of G, we have that {E) ^ C and = 

\n{Af^{E))\ ^ 01^1 with a constant 0 > 0 independent of ^ e G. For ^ e G, we 
define the mesh-size function 

G 3 X I—> hg{x) := |iY^({x})| ^ . 

For X G interior(£’), this coincides with the usual definition hg{x) = \E'^^'^ ='-hE- 
The mesh-size function is monotonic decreasing under refinement. 

We call the {d — 1 (-dimensional sub-simplexes of any simplex E e Q, whose 
interiors lie inside G, the sides of g and denote the set of all of them by S{g). For 
S G 5(0), we define hs := and observe for x G 5 that chs ^ hg{x) ^ Chs-, 

with constants 0, c > 0 depending solely on the shape-regularity of g. 

3.2. Finite element spaces. Denote by the space of polynomials of degree 
at most m G N. For a given grid 0 G G and certain subspaces Q c L'®(G) and 
V G 0(G), the finite element spaces are then given by 

V(0) := {y G V : V\e o g Py, F g 0 and Flan = o} , (3.1a) 

Q(e) := {q e Q : Q\e o g Pq, E g 0 } , (3.1b) 

where Py c W^’'^{E)‘^ and Pq c L^{E) are finite-dimensional subspaces Pi G 
Py G Pf and Pq G Pq G Pj for some £ ^ j G N. For convenience, we introduce the 
space of piecewise polynomials of degree at most m G N over 0 by 

Pm(0) := {i? : G ^ K : R\e e P^, Eg 0}. 

Note that Q(0) c L'^{gi) n Pj(0) and since V(0) c 0o(G) n P^(0) it follows that 
V(0) c iyQ’”(G). Additionally, we assume that the finite element spaces are nested. 
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i.e., if Qt is a refinement of Q, then 

V(g) cV(g*) 


and <Q{g) cz Q(C/*). 


(3.2) 


Each of the above spaces is supposed to have a finite and locally supported basis; 
e. g. for the discrete velocity space this means that for Q e G there exists an Ng e N 
such that 

V(e) = span{y?,...,ye^} 

and for each basis function Vf, i = 1,..., Ng, we have that if there exists an if e ^ 
with Vf # 0 on if, then supp V® c U^{E). We introduce the subspace of Vo(t/) of 
discretely divergence-free functions by 

Yo{g) :={f6V": J Qdiv Vdx = 0 for all Q 6 Q"} 


and we define 


Qo(0) := {geQ(g) :^Qdx = 0}. 


It will be assumed throughout the paper that all pairs of velocity-pressure finite 
element spaces considered possess the following properties. 

Assumption 7 (Projector Ufiv)- We assume that for each ^ e G there exists a 
linear projection operator —» V(C/) such that, for all s e (1, oo), 

• preserves divergence in Q(t/)*; i.e., for v e Wq’®(G)‘^ we have 
I gdivuda;= j gdivU^j^udcc for all Q e Q(0). 

Jq Ja 


^div locally defined; i.e., for any other partition e G we have 




(3.3) 


for all V e Wg^’*(f2)‘^ and all if e ^ with Af^{E) c 

Ujiv is locally W^’^-stable; i.e., there exists a ci > 0, independent of g, such 
that 


[ dx ^ Cl [ |u| 4- hg |Vu| dx 

Je Ju‘S(e) 


(3.4) 


for all V e Wg^’®(f2)'^ and all E e g. Here we have used the notation ■ dx := 
1 ^ • dx for the integral mean-value over a measurable set B c |ii| ^ 0. 


As in |DKS13I IBBDRl^ IDR07] . the local VE^’^-stability property (3.4| implies 
global W^’®-stability, i.e., for each s e [1, cx)], there exists a > 0, such that 




for all V e Wo’"(G)‘^. 


(3.5) 


Moreover, since V(^) contains piecewise affine functions, we have the following 
interpolation estimate. For all s e [1, oo], there exists a Cs > 0 such that 


/ 1^ - ^div’*^r + ^ “ '^^div'*^r da; < |u| 

J E 




(3.6) 


for all E e g and v e Wg^''’'^’®(G)'^, 0 ^ e ^ 1. 

As a consequence, we conclude the following result for weak limits in nested 
spaces. Before stating the result, we adopt the following notational convention: we 
shall write A < if to denote A ^ C ■ B with a constant C > 0 that is independent 
of the discretization parameter h. 
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Proposition 8. Let W^o’®(r2)‘^, s e (l,oo), be such that ^ 0 weakly 

in as k CO and let {Gk}kevi G be a sequence of nested partitions of 

LI, i.e., Gk ^ Gk+i for all fc e N. Then, 

^div^fc ^ 0 weakly in WQ^’®(ri)‘^ as k ^ co. 


Proof. Thanks to the uniform boundedness ( |3.5| ) of the sequence of linear operators 
- V(Cy,.) cz }fceN, we have that there exists a not rela¬ 
belled weakly converging subsequence of in Wg^’'^(n)‘^. By the compact 

embedding Wo^’'’(r2)‘^ L^{Ll)‘^ the sequence converges strongly in 

L^{Lty. Thanks to the uniqueness of the strong limit, it suffices to identify the limit 
of {‘^div'*^fe}fceN in L^{nY. To this end, we introduce the sets 

Gt ■■= n Gj and G: := {E e G+ : e t?+}, 


i.e., N^^{E) = J\f^>‘{E) for all j ^ k and E e G'^. For j > k, we consider the 
decomposition 

f^div^i ^ (‘^div'*^i)Xn(S + ) + (‘^div'*^i)Xn(efc\d,l')' 

For the latter term, we have according to ( |3.6| ) that 


lIVVo 


< 


^QkXu(g, 


k\Qt) 


llVt!, 


i'^3 f^div'*^j)xa(et,\e+) ^ ^ ^QjXu{gk\gf) 

Here we have used the monotonicity of the mesh-size under refinement in the last 
step. It follows from jMNSOSI Corollary 4.1 and (4.15)] that l^gkXQ(^g,,\g+) 

0 as fc ^ 00 . Thanks to the shape-regularity of G, this readily implies that 

= 0. (3.7) 


L°°(a) 


lim 
fe—>00 


^^GkXuig,\gf) 

By the compact embedding Wg^’*(n)'^ L®(r2)‘^ we have that Vj 0 strongly 

in L^{Ll)‘^ as j oo. Combining these observations, we conclude that for any e > 0 
there exists a iFf > 0 such that 


Pdiv^i)xa(efc\e^ 


^ e 


for all j ^ k ^ iFj. 


(3.8) 


We next investigate the term (3div'*^i)Xn(e+)- Thanks to the definition of G^ and 


(3.3) we have 


(fi'di'v^^j)ln(a+) = (^&v^j)la(e+) ^r all j > k. 


Since a linear operator between two normed linear spaces is norm-continuous if and 
only if it is weakly continuous (cf. Theorem 6.17 in |AB06| . for example,) we deduce 
that, for fixed fc e N, we have 

^ 0 weakly in W'^'‘‘{Ll{G^))‘‘ as j ^ go. 

By the compact embedding lFQ’*(f2)'^ L®(S7)‘^ this implies that 


(‘^div^i)Xn(e,l) 


strongly in L^(H)‘^ as j —» co. 


Together with (3.8), we have, for all j ^ k ^ K,,, that 

^'^div'^i)Xn(gk\gf) 


1 , . 

‘'div 


< 


^ e + 


(^div'*^i)Xn(g+) 


(‘^d(v'*^f)Xn(g + ) 

e as j 


00 . 


Since e > 0 was arbitrary, this proves the assertion. 


□ 


Next, we shall introduce a quasi-interpolation operator, which will be important 
for the treatment of the, generally non-polynomial, stress approximation. 
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Assumption 9. We assume that for each Q B G there exists a linear projection 
operator Ilg : such that Ilg is locally stable, i.e., there 

exists a c > 0, depending on such that 


[ ing^l dec ^ c /" I S'I da: for all S e 

Je Jus(e) 


lE JUS{E) 

As a consequence, this implies that 
lina^L^c. ||S||, 


for all S G L"(n) 


dx d 


(3.9) 


with a constant depending on Qq and s; compare also with (3.51. 

Assumption 10 (Projector 3 q). We assume that for each Q e G there exists a 
linear projection operator Uq : L^iTl) QiG) such that, for all s' e (l,oo), is 
locally stable, i.e., there exists a c > 0, independent of G, such that 


da; ^ 


Iu5(e) 


|g| da; for all q G and all E B G- 


We may argue similarly as for to deduce that 


II^Q^II^ < Cs llgll^ and / da;^Csh^^|g|^,,,(^g(£;)) (3.10) 

J E 

for all EbG and q B W'^’'^(n), 0 ^ e ^ 1. 

As a consequence of ( |2.3[ ) and Assumption [t] (compar e also with ( |3.5[ )) the fol¬ 
lowing discrete counterpart of ( |2.3[ ) holds; see |BBDR12] . 

Proposition 11 (Inf-sup stability). For all s, s' G (l,oo) with ^4-^ = 1, there 
exists a /?s > 0, independent of G ^ G, such that 
Q div y da: „ „ 

- ^PsWQl, for all QBQoiG). 

04VeV(a) IWIll.s 


Thanks to the above considerations, there exists a discrete Bogovskii operator, 
which has the following properties; compare also with |DKS13l Corollary 9]. 


Corollary 12 (Discrete Bogovskii operator). The linear operator S® := o 05 : 
divV(^) ^ V(^) satisfies 

Jf^HQdx 


div(»®il) = H and /3„ ^ ^ sup 


QeQiG) 


IIQL- 


for all E[ B divV(0) and s B (l,oo), with a positive constant /3s, independent of 
GeG. 

Moreover, let {t/fcjfceN G be a sequence of nested partitions ofGl, i.e., Gk+i ^ Gk 
for all fc G N, and let Vfc G Y{Gk) be such that ^ 0 weakly in iyQ’*(n)‘^ as 
k —>■ CO. ITe then have that 


divy^ ^ 0 weakly in Wg^’^(n)‘^ as fc —> oo. 

Proof. The claim follows as in |DKS12I Corollary 10] after replacing |DKS12I Propo¬ 
sition 7] by Proposition here in the proof. □ 

As in |Tem84| . we wish to ensure that the discrete counterpart of the convec¬ 
tion term inherits the skew-symmetry of the convection term. In particular, upon 
integration by parts, it follows that 

— (u (g) tu) ; Vhdx = / {vG)h) : Vw + (divu)(r(j • h) da; (3.11) 

J Gl J Q 

for all v,w,h B 'D(Gl)‘^. The last term vanishes provided that diva? = 0, and then 
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It can be easily seen that this is not generally true for finite element functions 
Fe even if 


f Q div F dx = 0 for all Q e Q(G), 

Jn 


(3.12) 


i. e., if F is discretely divergence-free. However, we observe from (3.11) that 
— / (v^w) :Vhdx = ^ / (v^h): Vw — (n (x) in) :Vhdx=: B[v, w, h] 

Jq Jq 

(3.13) 

for all v,w,he extend this definition to lF^’®(n)‘^ in the obvious 

way and deduce that 

B[v, n, n] = 0 for all v e 1F^’®(H)'^. (3.14) 

We further investigate this modified convection term for fixed r, r' e (l,oo) with 
i -I- ^ = 1; recall the definition of r from (2.11. We note that f > 1 is equivalent to 
the condition r > In this case we can define its dual f' e (1, cx)) by 4 -|- ^ = 1 
and we note that the Sobolev embedding 

holds. This is a crucial property in the continuous problem, which guarantees that 

f (i;(2)in) : V/ida; ^ c ||n||^ ^ ||in||^ ^ (3.16) 

for all v,w,he lF^’“(n)‘^; see |BGMSG09 . Because of the extension (3.13) of the 
convection term to functions that are not necessarily pointwise divergence-free, we 
have to adopt the following stronger condition in order to ensure that the trilinear 
form is bounded on x x (H)'^. In particular, let 

r > 2 ^, in order to ensure that there exists an s e (1, oo) such that y + ^ + j = 1- 
In other words, we have for v,w,h e lF^’”(n)'^ that 

/ (divn) (in ■h)dx^ l|divi;||^ ||in|| 2 j= \\h\\^ ^ c ||i;||i ||in||i \\h\\^ ~,, 

Jn 

with a constant c depending on r, and d. Here we have used the embeddings 
(3.15) and (TlY ^ L®(r2)‘^. Consequently, together with (3.16) we thus obtain 


B[v, in, ft.] ^ c ||n||^_^ ||in||i_^ ||ft||^_-, . 


(3.17) 


In view of (3.13), for v = (ni,... , 1 )^)" e Wg’’^(n)“, the convective term can be 
reformulated as 


B[n, n] • in da; = S[n, n, in], in e IFq’’^( fl)'^ 


(3.18) 


where B[n,i;] e is defined by {B[v,v])j = 

j = 1,..., d. In particular, for n = F e V(^), we have that B[F, F] e V 2 e.-i{G)‘^- 

Example 13. The following velocity-pressure pairs of finite elements satisfy As¬ 
sumptions 0 and [To] for d = 2,3 (see, e.g., [BBDRl^ IGLOll lGS03p : 

• Lowest order Taylor-Hood element; 

• Spaces of continuous piecewise quadratic elements for the velocity and piecewise 
constants for the pressure ([BF9T] §VI Example 3.6]). 

We note that the MINI element and the conforming Crouzeix-Raviart Stokes ele¬ 
ment do not satisfy the nestedness hypothesis stated in (3.2). 
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Remark 14. The boundedness of the trilinear form B[-, •, •] stated in (3.171 re¬ 
quires that r > In |DKS12| and |DKS13| the set of admissible values of r 

was the same range, r e as in the existence theorem for the continuous 

problem in |BGMSG09| ; however, for r e ths finite element space for the 

velocity was assumed in |DKS12| and |DKS13| to consist of pointwise divergence- 
free functions, whose construction is more complicated. For simplicity, we shall 
therefore confine ourselves here to the limited range of r > so as to be able to 
admit standard discretely divergence-free (cf. (3.12[)^ finite element velocity spaces. 


3.3. The Galerkin approximation. We are now ready to state the discrete prob¬ 


lem. Let {V(0), Q(^)}ggG be the finite element spaces of Section 3.2 

For n G N and Q e G we call a triple of functions {Ug, Pg) e V(C/) x Qo(S) a 
Galerkin approximation of (2.6 1 if it satisfies 


f -.DV+ B[U2,U2]-V - P^dwVdx = [ f-Vdx, 

Jq Jq 

[ QdivtJg da: = 0, 

Jq 


(3.19) 


for all V G V(^) and Q G Q(^;). 

Restricting the test-functions to Vo(5) the discrete problem (3.231 reduces to 
finding Ug G Vo(0) such that 

[ S^{-,DU^) ■.DVdx + B[U'^,U^,V]= f f-Vdx (3.20) 

y f2 y f2 


for all V G Yo{G). Thanks to (3.141, it follows from Assumption and Korn’s 
inequality (2.5) that the nonlinear operator defined on Vo(^) by the left-hand side 


of (3.20) is coercive and continuous on Vo(0)- Since the dimension ofY{G) is finite, 
Brouwer’s fixed point theorem ensures the existence of a solution to (3.20). The 


existence of a solution triple to (3.23) then follows by the discrete inf-sup stability. 
Proposition [m Of course, because of the weak assumptions in the definition of 


the maximal monotone r-graph, (3.23) does not define the Galerkin approximation 


{Ug, Pg) uniquely. However, supposing the axiom of choice, for each n G N, ^ G G, 
we may choose an arbitrary one among possibly infinitely many solution triples and 
thus obtain 


{{U^,P^,S-{;DU^))} 


neN.geG’ 


(3.21) 


From (3.23) we see that Ug is discretely divergence-free and thus, thanks to (3.20) 
and (3.14), we have that 

[ S^{;DU^):DU^dx = (f, i|/||_i,, WU'^l,. 

Jq 

The coercivity of 5” (Assumption]^ and Korn’s inequality (2.5) imply that the 
sequence is bounded in the norm of Wg^’’^(n)'^, independently of ^ G G and 

n G N. This in turn implies, again by Assumption the uniform boundedness of 
S^{-,DUg) in L’’ {Il)‘^^‘^. In other words, there exists a constant Cf > 0 depending 
on the data /, such that 


\\U2l^+\\S-{;BU2)l,^Cf, 


for all ^ G G and n G N. 


(3.22) 


For the sake of simplicity of the presentation, if there is no risk of confusion, we 
will denote in what follows S^{I)Ug) = jS”(-,DC/g). 
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Remark 15. An alternative formulation of (3.201 is follows: find a triple {Ug, Pg, Sg^ e 
y(Q) X Qo(S) X sueh that 


Sg :Dy + B[J7g,i7g]divFdx= / f-Vdx, 

J fi 

[ QdivU^ dx = 0, 

in 


(3.23) 


/ S'^:Ddx= / : Ddx, 

I n in 

\dxd 


for all V e V(^), Q e Q(^), and D e In particular, if we define 

iig: by 


dxd 


/ rigs':Z3da;= / S:Ddx, for all D e V(-i{Q) 
in in 

then Assumption^ can be easily verified and Sg may take the role of AgS'^fDUg) 
in the subsequent analysis. 


3.4. Discrete Lipschitz truncation. In this section we shall recall a discrete 
counterpart of Lipschitz truncation, which acts on finite element spaces. This dis¬ 
crete Lipschitz truncation is a composition of a continuous Lipschitz truncation 
with a projection onto the finite element space. The continuous Lipschitz trunca¬ 
tion used here is based on results from jDMSOSI IBDF121 IBDS13| . which provides 
finer estimates than the original Lipschitz truncation technique proposed by Acerbi 
and Fusco in |AF88| : for details consider |DKS13| . 

We summarize the properties of the discrete Lipschitz truncation in the following 
result. Similar results for Sobolev functions can be found in |DMS08| and |BDF12| . 


Proposition 16. Let 1 < s < oo and let {Ek}keN a sequence such that for all 
k e N we have e V(^fc) for some Qk G G. In addition, assume that {Ek}kefi ^ 

converges to zero weakly in Wq'^, as fc —» oo. 

Then, there exists a sequence {Xk,j}k,jeN ^ K with 2^^ ^ Xkj ^ 2^^^ and 
Lipschitz truncated functions Ekj = Ek^x,.., k,j e N, with the following properties: 

(a) Ekj eY{gk); 

(b) < c forl<s^ oo; 

(c) II ||oo ^ C XkJ i 

(d) Ekj —> 0 in L’^{n)‘^ as k ^ oo; 

(e) YEk,j —* 0 in L’^{ny^‘^ as k ^ oo; 


(f) For all k,jeN we have ||Afej- X{ 


Ek¥^Ek,j} 


^ c2- 




The constants c appearing in the inequalities (b), (c)| and (f) depend on d, LI, Py 
and the shape-regularity of {Qk}keN- The constants in|(b)| and|(f)| afeo depend on s. 


Proof. The proof is exactly the same as that of |DKS12l Theorem 17 and Corollary 
18], replacing [DKS121 Proposition 7] by Proposition]^ □ 


4. Error Analysis 

4.1. Graph approximation error. In order to quantify the error committed in 
the approximation of the graph A{x), x B LI, we introduce the following indicator. 
For D e L^{Liy^‘^, S e L^'{Liy^‘^, we define 

£a{D, S) := [ inf |D - + 1^ - dx. (4.1) 

Jn (S,rT)€A(x) 

The following result shows that this indicator is well-defined. 
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Proposition 17. Let D 6 and S B L'' then, the mapping 

^ “ ‘^ 1 ’' + “ '^ 1 '' 

(d,<T)ev4(a;) 

is integrahle. Moreover, there exist D e and S B L'^ such that 

{jD{x), S{x)) B A{x) for a.e. x B LI and 


£a{D,S) 


D-D 


r 

+ 


s-s 


da;. 


Proof. The first claim is an immediate consequence of the second one. The second 
assertion follows from |AF09l Theorem 8.2.11] by observing that the mapping 


11 X X 9 (x; (.5, <t)) - |Il(x) - <5r + \S{x) - cr^' 

is Caratheodory, i.e., x |Z3 + |S^(x) — o'l’' is measurable for all cr e 

and \D{x) — J]’’ + |S'(x) — cr]’’ is continuous for a.e. x B Lt. □ 


4.2. A posteriori finite element error estimates. In this section we shall prove 
bounds on the residual 

n{U^,P^,S^{DU'^)) = (7^P''"(^7g,P^,S'”(D^7g)),7^''^(^7g)) e W-^’^{L}fxLl{Ll) 

of ( |2.6| ). In particular, for {v,q,T) B lTQ’^(n)‘^ x Tq(0) x we have 

(n{v, q, T), (w, o)} := (nP^%v, q, T), w) + (U'^iv), o) 

:= / T : Dm + B[u,u] • to — qdivtu — / • Id — odivudx, 
an 

(4.2) 

where (m,o) e x L’' (n)/]R. Although for the sake of simplicity we re¬ 

strict ourselves here to residual-based estimates, we note that in principle other a 
posteriori techniques, such as hierarchical estimates, flux-equilibration or estimates 
based on local problems, can be used as well; compare with |MSVn7l ISieOQj . For 
n e N and t/ G G let {Ug, Pg) b Y{Q) x Qo(fJ) be the Galerkin approximation 
defined in (3.231. We begin with some preliminary observations. 

The first part of the residual in (|4.2|), q, T) B provides infor¬ 


mation about how well the functions v,q,T satisfy the first equation in (2.6). For 
the second part, we have TZ"^{v) B (Lg(n))*. We note that the space {Lq(LI))* is 
isometrically isomorphic to L’' (D)/K, which is, in turn, isomorphic to Lq (D) since 
r G (l,oo). The term (TZ'‘^{v), o') provides information about the compressibility 
of V. 

We emphasize that Tl{v, q,T) = 0 if and only if q,T) = 0 and TZ"^{v) = 0, 

but that a vanishing residual itself does not guarantee that (Du(x),T(x)) G A{x) 
for almost every x bLI. For this, additionally £’_4(Dr;,r) = 0 is needed. 

For the rest of the paper let t and i be such that 


2 d 

d + 1 


< t < r 


t = r 


and 

and 


1 dt 


2 d-t’ 

i = t' = f = r', 

3d 


if r ^ 


3d 


d4-2’ 
otherwise. 


(4.3a) 

(4.3b) 


Note that (4.3a) implies that if r ^ then t < r and t 


< r. 


Lemma 18. The triple S') G x Lq(S1) x L'^ (fl) 

(2.2) if and only if 

Tl{u,p, S) = 0 in X L* (D) 


dxd 


is a solution of 


and S') = 0. 
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w^/{nY 


X L‘'(ri)/]R is dense in {^Y ^ 


Proof. Thanks to the fact, that 
L’’(r2)/K we have that Tl{u,p,S) = 0 in W~^'*'{PtY x Lq{PI) is equivalent to 
TZ{u,p,S) = 0 in {YY x L^ip). This is, in turn, equivalent to the fact 

that the triple (m,p, S) satisfies the system of partial differential equations (2.2). 

On the other hand we have that (Dtt(a;), S{x)) e A(x) for almost every a; e O if 
and only if £pDu, S) = 0, and that completes the proof. □ 


Note that Lemma does not provide a quantitative relation between the error 
and the residual. Even for simple r-Laplacian type problems, such a relation requires 
complicated techniques and problem-adapted error notions (e.g. a suitable quasi¬ 
norm); cf. |DK08I IBDK12I rLYfll| . However, because of the possible nonuniqueness 
of solutions to (2.21, such a relation cannot be guaranteed in our situation. We shall 
therefore restrict the a posteriori analysis to bounding the residual of the problem 
instead of bounding the error. 

Recalling the quasi-interpolation Hg from Assumption as well as the represen¬ 
tation of the discrete convective term in ( |3.18[ ), we define the local indicators on 
E e Q as follows: 


t 

i,dE 


£f^{Ug,PY S^pUgfE) := \\hg{- dlYUgS'^pUg) + B[Ug,Ug] + VPg 

+ pgS'"(Dt7g)-Pg"idl 

+ ||S"(DC/e)-ng5"(Di7g)||;-, 

4^(C/e;P) := ||divJ7g||;^, 



(4.4a) 

(4.4b) 


and 


£g{U^,P^,S^{DU^y,E) :=£^g^YUg,Pg,S^\DU^yE)+£^YU^-E). (4.4c) 

Here, for S e S{G), [[•]] |s denotes the normal jump across S and [[•]] := 0. 

Moreover, we define the error bounds to be the sums of the local indicators, i.e., for 
Ai cz Q, we have 


£g{U^;M) 


= 2 £:P'"(i7g,P^\5”(Dt7g);P), 

EeM 


div 


t 

t-.niM) 


and 


£g{U^,P§,S^{BU’^)) 


:=£^g^YUG,Pg\SYDU-g))+£'^YUg) 

: = (t/g, Pg", 5" (D C/g); g) -f 4^ ( C/g; 0). 


Theorem 19 (Upper bound on the residual). Let n e N and Q e G, and denote by 
([/”, Pg) e Y{G) X Qo(0) o, Galerkin approximation of (3.23). We then have the 
following bounds: 


||7^P^"(C/g,Pg^^"(D^7g))|| 




o€L*'(f2)/K 


sup (P'4«7g),T 


infceB ||o- c|L, 


<Ci£f^(i7g,Pg",5"(Dt7g)) 

= 4(c/g)'^‘. 


i/t 


(4.5a) 

(4.5b) 


The constant Ci > 0 depends only on the shape-regularity of G, t, and on the 
dimension d. 
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Proof. The assertions are proved using standard techniques; compare e.g. with 
[Ver96lfAOOO| . For the reader’s convenience we sketch the arguments. For arbitrary 
{v,q) e Wq’^ X L* (ri)/]R with = \p\^, = 1 we deduce from (3.23) that 

= [ ngS"(DC/g) : D(u - + B[J7g, C/g] • {v - - f ■ (v - dx 

Jn 

- f PSdiviv-J%^v)dx+ f (S'"(Dt/g)-ne5"(Di7g)) :D(u-ae^i;)dx. 
Jn Jn 

Thanks to ( |3.18 1, local integration by parts and using Holder’s inequality, we obtain 
<7^P^^(^7g,Pg^5"(D^7g)),u> 

2 { I- divng5”(D[/g) + B[i7g, C/g] + VPe" - f^E ll« - 


EeQ 


< 


+ 2 \\lUoS-{BU-g) - PSmi.oE ||« - 
+ ||S"(DC/g) - ne5"(Di7g)||-^ l|Du||p_^ } 

E { ll^e( - divneS"(DJ7g) + B[c/g,c/g] + VPg" - f)\\l^ 


EeQ 


lit 

\\i,eE 


+ \\hy* pgS'"(D[/g)-Pe"idl| 

+ ||S"(DC/g) - ne5"(D«7g)||;;^ hlliy,. 


Here, in the last inequality, we used the stability of (see (3.4|), a scaled trace 
theorem, and the interpolation estimate for in (3.6), as well as the finite over¬ 
lapping of patches and a scaled trace theorem. 

To prove the bound (4.5b), we first conclude from J^ldivC/gda; = 0 and 
Holder’s inequality that for all c e K, we have 

f odivUgdx= f (o — c) div Lfg da; ^ ||div Lfgl ||o — c||(, . 

Jn Jn 

Taking the infimum over all c e K and then the supremum over all o e L’' (f2) proves 
in (|4.5b[). In order to prove we observe that 


£'g=(J7g) = [ divJ7g |divi7g|‘-2divi7gda; 

Jn 

^ sup /p-(i7g), 

o€L*'(n)/R \ infceR ||o c||(, 


div i7g|‘-^div Lrg| 


Together with the definition of Sg(Ug) and noting that 

|||divC/g|‘-2divi7g||^, = ||divC/g||;-' =4‘=(t7g)'-^ =£:g^(i7g)- 


this yields (4.5b). 

Corollary 20. Under the conditions of Theorem ]! 9[ we have 

(np^%u2,p^,s^fDU2 )),v) 

^ Cl ^ £-P^^(Cg,Pe",5”(DCg);P)'/‘||Vu|| 

i',ue{E) 

EeQ 


□ 
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and 


Eeg 

for all V e Wq'^ {ilY qB L* (fi). 

Theorem 21 (Lower bound on the residual). Under the conditions of Theorem^H 
we have 


C2£-P^^(«7g,Pe",5"(D«7g)) 


1/i 


< 




i/t 


(4.6) 


The constant C 2 > 0 solely depends on the shape-regularity of Q, t, and on the 
dimension d. The oscillation term is defined by 


osce(t7g,S"(DC/g)) := 2 osc(C/e,5"(Di7g),P) 

EeO 


Eeg 


min 

f 


ho{f-fE)\\lE 


+ ||5"(Di7g)-ne5"(D[/g)||;;^. 


Proof. Let E e Q and let S e S(C/), i.e., there exist Ei,E 2 e Q, Ei E 2 , such that 
S' = Pi n i? 2 - Let G P 2 £-i arbitrary; for convenience we use the notation 


Re :=- divneS”(D[/g) + B[U'^,U^] + - /^ G P^,_i 


and 


Js := lngS^{BUl) - Pg" id] Is e P^^^ where m = max{£ - 1, j}. 


It is well known that there exist local bubble functions Be, £ M^o’”(r2), such 
that 

0 ^ 6 ej Bg < 1, supp Be = E and supp Bg = ojg := Ei u ^ 2 - (4.7a) 

Moreover, we have that there exist pE e P 2 ^_i and pg G with ||pe||(/ e ~ ^ ~ 

llpsllt' s> S'^ch that 


/ ReBePE^X, ||V(6£;P_e)||p ^ ^ C ||/lg^p£;||-, , 

’ JE 


Il'^sllt.S jg '^sl’SPS^X, l|V(6sPs)|lt',cc;s ^ ^ 
and ||^SPs|lt/,^s < C h]l^ pg _ ^ ; 


L-i/f' 

«g PS 


i’,S ' 


(4.7b) 


compare e.g. with |Verl3l Chapter 3.6]. Here the constants only depend on r, the 
polynomial degree oi Re, respectively Jg, on the shape-regularity of Q, and on the 
dimension d. Hence, for the element residual, we conclude that 

WReWie < C'{<pP^^(«7g,Pe”,S”(Dt7g)), BePe) 

+ <ngS"(D[/g) - S"(Dt70), D{BePe)) + 11/ - fEWi.E }. 
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where we have used Holder’s inequality and that 0 ^ ^ 1. Together with a 


triangle inequality and (4.7bI, this implies that 


\\hg{- divngS"(DC/g) + B[i7g,i7g] + VPg" - /)||-^ 

< c{<pP^nC/e,Pe",5"(DC/g)), hgSEPE) (4.8) 

+ \\UgS-{DU-g) - S”(DC/g)||-^ -f \\hg{f - fE)\kE }• 


For the jump residual, we deduce from (4.7al and integration by parts, that 


UsWis^C f pg 5" (D iJg) - Pg” id] Bspsds 

’ JS 

= c{<7^P^^(^7g,Pg^S”(D^7e)), Bsps) 

+ y f (divng5"(Dt7g)-B[C/e,J7e]-VPg" + /)6sPsdx}. 


Therefore, we obtain, with (4.7b), Holder’s inequality, and (4.8), that 


/if l5"(Dt7g)-Pg"id] 


i,S 


^ c{(pP^^(C/a,Pg",5"(DC/e)), hfssps] 

+ E [<PP'^(C/g,Pg",S'"(D[/e)), /ig^pB.) 

2 = 1,2 ^ 

+ ||ng5"(D17g) - S"(DC/^)||-^^ -f \\hg{f - fEjllE.]}- 


(4.9) 


We define := ||/ig(-divng5"(D«7g) + B[J7g, [/()]+VP^-/) ||‘-f E e G, 


\\t-l 


and Ps ■■= pf [[ngS'"(DJ7g)-P^id]]|_^, S 6 S{G). Then, combining (|0) 


and (4.9) and summing over all P e 5, 5” e 5(^), yields 


ffpt7g,Pe",S'"(Dt7g)) 


= ^ ai5pg(-divngP"(DJ7g)+B[[/g,[/g] + VPg"-/)||-^ 


Eeg 


seSig) 


y Ps /if lP”(DJ7g)-Pg"id] 


t,S 


< 


c\ /pP^^(i7g,Pg",S'"(D[/g)), y(aE+ y Ps)hg6EPE) 

(. \ Eeg ^ SadEnQ / 

+ /pP^^([/g,Pg",P"(DC/a)), 2 Pshfbsp^ 

\ SeS{g) / 


yy^EE y Ps)of^c{iJi,s-{^ui),E) 

EeG SddEnQ 


1/i 


Here we have used that f ^ e E e G are arbitrary in the last step. 
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Thanks to the fact that the supp 6^, E e Q, are mutually disjoint up to a nullset 
together with (4.7bI, we have that 


SczdEnQ 


Xj + X Ps^^ihgbEPE) 

= X (“^ + X [ \hG\^{6EPE)\* dec 

ZT TP^n J E 


E SezdEnO. 

^c[Y,Je+ X 4') 

Se5(e) 


E 

rpde 




1/i' 


where we have used that each element E has at most (d + 1) sides S e S{Q) with 
S cz dE. The constants C depend only on the shape-regularity of Q. Analogously, 
we conclude from the fact that only finitely many of the supp Bg, S e S{Q), overlap, 
that 

( X f^shf^sps)\[^c Y, /3*'/ Ihfy^psf Y, P'' 


SeS(a) 


S€5(e) 


S^S(G) 


,1/t' 


Combining Holder’s inequality with similar arguments yields for the last term 


X(«s+ X /35)osc(l7g,S"(DC/^);A) 




S^zdEnQ 


1/i 


2 4')'^Xsc(C/g,5"(D[7g))i/‘ 


EeG 


SeSiG) 


^ C£f\Ul,P^, 5”(D«7g))^/‘' osc(«7g, 
Altogether, we have thus proved that 

+ osc(«7e, 5"(D17g)) 

This is the desired bound. 


□ 


The following result states the local stability of the error bound and is referred 
to as local lower bound in the context of linear elliptic problems. 


Corollary 22 (Local stability). Suppose the conditions of Theorem 19 and let 
A4 c Q; then, there exists a constant C depending solely on the shape-regularity of 
Q, t, d and El, such that 


i/t 


( \nP^%Ul,P^, + osc(C/g, S^{BUl)-M) 

^ \Pg\\i,.p>./m) + \Pg\\,,uo/M) + PokuHM) + WfkuHM) + k 


< C 


1 /i 


r,U'S(M) 


)■ 


Proof. The first bound follows as in the proof of Theorem In order to prove 
the second bound, let v e VLq’* {U^{M)y. We then have with Holder’s inequality 
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and ( |3.17[ ), with t and i instead of r and r, that 
<7^P^^(C/g,Pe^5"(D^7g)), i;> 

= [ UgS^iDU'^) ■.Dv + B[U'^,U^]-v-P^divv-f-vdx 

Jq 

< c( ||n5S”(D£/5)||, 


P. 


s \\i-Me(M) 


+ II /lit; 


i-,US{M) 


lll,t';W6(A1) • 


Note that in the case of (4.3aI, we have 

1 dt 


t 


~ 1 dt 1 dt t .t 

t = -- < -- = r- ^ r - = -- =: s'. 

2 d — t 2 d — r r r r — 1 

Hence, Holder’s inequality, the stability of Bg (Assumption and Assumption 
yield 


H 








Bu: 


n II I 




r'M9{M) ■ 

The oscillation term can be bounded above similarly, and the assertions follows. □ 

Remark 23. Corollary \2t^ states the stability properties of the estimator, which 
are required in order to apply the convergence theory in |Sie09l lMNS08| .' compare 
with |Sie09l (2.10b)], for example. The stability of the estimator is also of im¬ 
portance for the efficiency of the estimator. If Corollary fails to hold, it may 
happen that the a posteriori error estimator is unbounded even though the sequence 
of discrete solutions is convergent; in particular, divS'"(D{7g) need not belong to 
U" (E) when 1 < r < 2. This problem already appears in the a posteriori analy¬ 
sis of quadratic finite element approximations of the r-Laplacian, or the r-Stokes 
problem (cf. |BS08) ) for 1 < r < 2. In order to avoid this, we use ngS"(DC/g) 
instead of S^fDUg) in the element residual (4.4a). This is compensated by the 

term ||5'"(D[7g) — ngS”(Di7g)||~ in the a posteriori bounds (4.Sal and (4.6); cf. 
Appendix\^for further details. 


5. Convergent Adaptive Finite Elements 

This section is concerned with the proof of convergence of an adaptive finite 
element algorithm for the implicit constitutive model under consideration. 

5.1. The adaptive finite element method (AFEM). In this section, we shall 
introduce an adaptive finite element method for ( |2.2| ) . The details of the subroutines 
used in the process are listed below: 

The routine SOLVE. We assume that for arbitrary n e N, Q e G, the routine 
S0LVE(n, Q) = {Ug,Pn, S'^{-,DUg)) computes an exact solution (Ug,Pn) e V(^) x 
Q(^l) of 

The routine MARK. For a fixed function g : Kj —> Kj, which is continuous at 0 with 

g(0) = 0, we assume that the set At = MARK(f£g (Ug, Pg, S"'(DUg); 

satisfies 

max{£g(U^,P^,S"(;DU^);E) : E e g\M} 

^g(max{£g(U^,P^,S^(;DU^);E) :EeM}). 


(5.1) 
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Algorithm 1 AFEM 

Let fc = 0, no = 1) and let Go be a given partition of fl. 


1: 

loop forever 


2: 

3: 

let Sk = S'”". 

(Uk,Pk,Sk(-,DUk)) = S0LVE(nk,Gk) 


4: 

compute {£g,^(Uk,Pk,Sk(-,DUk);E)}^^g^, and £:_^(DUk, Sk(-,DUk)) 

5: 

if £g,(Uk,Pk,Sk(-,DUk)} 5^ £A(DUk,Sk(;DUk)) then 

6: 

Mk = MARK({£;g,^(Uk,Pk,Sk(-,DUk);E)}j^^g^ 

,Gk) 

7: 

Gk+i = REFINE(Gk,Mk) 

% mesh-refinement 

8: 

'^k+1 ~ ’^k 


9: 

else 


10: 

Uk+l = Uk + f 

% graph-refinement 

11: 

end if 


12: 

k = k + 1 


13: 

end loop forever 



Hence the marking criterion guarantees that all indicators in Q are controlled by 
the maximal indicator in Ai. Note that this criterion covers most commonly used 
marking strategies with g{s) = s; cf. |Sie09l IMNS08] . 

For the definition of the routine REFINE see Section IXTI 

For the sake of simplicity of the presentation, in the following, we will suppress 
the dependence on x in our notation and write Sk(DUk) = Ski’y^Uk) if there is 
no risk of confusion. 

5.2. Convergence of the AFEM. Let {Gk}kefi cz G be the sequence of meshes 
produced by AFEM. For s e (l,oo], we define 

:= U and := c Lg(L!). (5.2) 

k^O k^O 

Lemma 24. Let {{Uk, Pk. SkiBUk ))}c x L5(L!) x be the 

sequence produced by AFEM; then, at least for a not relabelled subsequence, we have 

Uk Moo weakly in 

Pk Poo weakly in Lq (LI ), 

Sk{DUk) — Sao weakly in 

for some (moojPooj S'oo) £ x x Lg (fl). Moreover, we have that 
niUk, Pk, SkiBUk)) 7^(Moo,Poo, 5oo) weakly* in 

and 

<7^(Moo,Poo,S'oo), {v,q)} = 0 for all gG G V^. 

Proof. The proof is postponed to Section |6.1| □ 

Corollary 25. Let {{Uk,Pk, Sk{P>Uk))}^^^ c x Ll{Ll) x u'be a 

not relabelled subsequence with weak limit (mooiPooj S/ao) s x x L’’ as 

in Lemma Then, 


£Gk{Uk,Pk,Sk{P>Uk)) asfc^oo; 
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implies that 

7^(^lao,Poo, -Soo) = 06 


Proof. The upper bound, Theorem [I^ together with £gi^{Uk, Pk, Sk(DUk)) 0 
as k ^ CO, implies that 

n{Uk,Pk,Sk{BUk)) ^ 0 strongly in 

Thus the assertion follows from Lemmaand the uniqueness of the limit. □ 

Lemma 26. Let { (iJfe, Sfc(Di7fc)) c x L5(L!) x U'he a 

not relabelled subsequence with weak limit (moo,Pooj S^oo) S x x as 

in Lemma Assume that 

£A{P>Uk,Sk(DUk)) ^ 0 as k ^ CO; 


then, 


(DMoo(a;), 5'oo(a^)) s A{x) for almost every x e LI. 


Proof. The proof of this lemma is postponed to Section |6.2| below. □ 

Lemma 27. Assume that the sequence {nk}keN satisfies Uk —>■ N < co as k —>■ co. 
We then have that 


£g,^{Uk,Pk,Sk(DUk))0 as k^ CO. 


Proof. The proof of this lemma is postponed to Section [63] below. 


□ 


We further assume that the graph approximation is uniform with respect to the 
graph approximation indicator. 


Assumption 28. For every e > 0, there exists an A = N{e) 6 N, such that 

£:_ 4 (Di;,S'”(-,Dn)) < e for all v 6 and n> N. 

We note that this and Assumption are the only strong assumptions among the 
ones we have made; Assumption |28| is, however, only used in the proof of the next 
theorem, and is not required for any of the preceding results. 


Theorem 29. Let {{Uk, Pk, SkfDUk))} be the sequence of function triples produced 
by the AFEM. We then have that 

£A{DUk,Sk(DUk)) ^ 0 as k ^ CO 


and, for a not relabelled subsequence, we have that 


£gkiUk,Pk,Sk(DUk)) ^ 0 as k^ CO. 


Proof. We argue by contradiction. First assume that there exists an e > 0 such 
that, for some subsequence, we have that 


£A{^Uk„Sk,{^UkJ)>e, foralHeN. 

Consequently, by Assumption|^ we have that Uk^ = N, for some £o, A 6 N, and all 

there exists a not relabelled subsequence 


> 


Hq. Moreover, thanks to Lemma 24 
{{Uke, Pke, SkifDUkt))} that converges weakly in WS{n)‘^ x L'^{L1) x L^'(Lt) 
Combining these facts, we conclude with Lemma [27| that 

^GkAUke, Pki, Skt{P>Uke)) -♦ 0 . 


dxd 


In particular, there exists an £ > £o, such that £g,^^{Uki, Pkt, < e. 

Therefore, by line of AFEM we have that Uki+i = N + 1, a, contradiction. 
Consequently, we have (for the full sequence) that 


£A{OUk,Sk{DUk)) ^ 0 as /cGO. 


This proves the first claim. 
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Assume now that there exists an e > 0 such that we have that 
8g^{Uk,Pk,Sk{J^Uk))>e for all A: sN. 


(5.3) 


By the above considerations, there exists a /cq e N, such that £yi(DLffc, Sk{P>Uk)) < 
e for all fc > fco- Therefore, according to line of AFEM, we have that Uk = rzfeo 
for all k ^ kg. Consequently, Lemma 27 contradicts (5.3). 

Combining the two cases proves the assertion. □ 

Corollary 30. Let {(Uk, Pk, Sk(DUk))} be the sequence of function triples pro¬ 
duced by the AFEM. Then, there exists a not relabelled subsequence with weak limit 
{uaojPoo, Sao) G [LlY X Lq{LI) X L’’ such that 

8 A{^Uk,Sk{DUk)) ^ 0 and £g,^{Uk, Pk, Sk{DUk)) ^ 0, 

as k CO and (moojPoo solves (|2.2[). 


Proof. The claim follows from Theorem]^ Lemma|^ Corollaryj^ and Lemma[T8| 

□ 


Remark 31. ITe emphasize that even in the case when the exact solution of (2.2) 


is unique, we do not have that the statement of Corollary \3C\ is true for the full 
sequence. This is due to the fact that the finite element error estimator is not 
necessarily decreasing with respect to the refinement of the graph approximation. 
However, when the exact solution is unique, it is easy to select a converging sub¬ 
sequence with the help of the estimators; choose, for example, a subsequence, such 
that £g,,{Uki,Pki^ Ski{T>Ukf)) is monotonic decreasing in £. 

6. The Proofs of the Auxiliary Results 


6.1. Proof of Lemma 24 We recall (3.22) and observe that the spaces WQ’^{n)'^ 
and L’" r e (l,oo), are reflexive. Therefore, there exist ttoo G and 

So 3 s L’’ such that for a not relabelled subsequence we have 


Uk 


Ur, 


weakly in WQ’^(r2) 


( 6 . 1 ) 


and 


Sk{T>Uk) — Sao weakly in (6.2) 

as fc —» cx). The function Mqo is discretely divergence-free with respect to Q^, i.e.. 


q div Mqo da: = lim / q) div Ukdx = 0 for all g e ' 


fc —>-00 


This follows from (3.10) as in the proof of |DKS131 Lemma 19], replacing |DKS13I 
(3.5)] with the density of the union of the discrete pressure spaces in Q^. 
Moreover, using compact embeddings of Sobolev spaces, we have that 

^sG(l,f^), ifr<d, 

G (l,oo), otherwise. 

This implies, for arbitrary v e Wq’“( fl)'^, that 
B[Uk, Uk, u] ^ B[uao, Mqo, •u], 


Uk 


strongly in for all 


(6.3) 


or equivalently, 

B[Uk,Uk] - B[u^,u^] weakly in 


as A: —> cx); compare also with |DKS13l Lemma 19]. 
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We now prove convergence of the pressure. Thanks to (3.171, we have 

f PkdivVdx= f Sk{BUk) -.ddY + B[Uk,Uk] ■ V - f -V dx 
Jq Jn 

^ ||5fc(Dt7fc)t, iiDyt + c WUkWlr ||y||i^, 

for all V e Y{Qk)- By (3.22) and the discrete inf-sup condition stated in Proposi¬ 
tion]^ it follows that the sequence {Pk]k€N is bounded in the reflexive Banach space 
Lg(n). Hence, there exists a Poo e c Lq{^) such that, for a (not relabelled) 
subsequence, 

Pk Po 3 weakly in Lo(H). 

On the other hand we deduce for an arbitrary v e V® c; that 


Poo div V da; 


Pk div vdx = / Pk div jfPv da: 


= [ SkiBUk) : - f ■ + B[Uk, Uk] ■ da; 

an 


^ / Soo : Du da; + B[moo, Moo] • l: — / • da; 
an 

as A; ^ 00 , where we have used (6.2), the properties of together with the 
density of the union of the discrete velocity spaces in and the boundedness of 
the sequence {P/c}neN in Dhe assertion for all v e then follows from the 

density of V® in V^. 


6.2. Proof of Lemma 

bk G and S 

X e Q and 


□ 

According to Proposition for A; e N there exist 
k U such that {bk{x), Sk{x)) G A{x) for a.e. 


Dt/fe — Di 


-f 


5fc(DC/fc)-Sfe =£A{BUk,Sk{BUk))^d ask^ao. 

r' 

(6.4) 


Thanks to (3.221, the sequences {bk}keN and {Sk}keN are bounded in and 

L’’ (n)'^^'^ respectively. Since both spaces are reflexive, together with the uniqueness 
of the limits, we obtain that 

bk Du^ weakly in (6.5a) 

Sk — Sao weakly in (6.5b) 

as k ^ CO. 

Let S* : X ^ ®sym a measurable selection with (d, S*{x,S)) 6 A{x) 
for a.e. x e and thus (DMoo(a;), S*(a;, DMoo(a;))) G A{x) for a.e. x e Q ;compare 
with Remark]^ Consequently, for every bounded sequence {^fcjfceN £ L®(n) of 
nonnegative functions, we have (recall [(A3) [) that 


0 ^ limsup 
fc —>00 ^ 


Sk - S'*(Dmoo)) : (Dk - Duoo) (j^kdx 


= limsup / [Sk - S*{Duao)) : {bk - DuaA) 4>k dx 


fc —>00 


( 6 . 6 ) 


= limsup f {Sk{BUk) - S*{Duao)) ■■ {^Uk - Duoo) (pkdx, 

fc—>CX) Jq ---v--- 


= -.ak(x) 


where we have used (6.4) in the last step. We assume for the moment that ak 0 
in measure and therefore 
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for at least a subsequence of a^- Since is bounded in L^(n), we obtain with the 
biting Lemma (Lemma and Vitali’s theorem, that there exists a nonincreasing 
sequence of measurable subsets Ej c with \Ej\ 0 as j ^ oo, such that for all 
j e N, we have 


Ofc 


0 strongly in as fc ^ oo. 


This, together with (6.6l and (6.5), implies for all nonnegative 4> ^ 
L’^iyi) (extend (p by zero on Ej) and each fixed j e N, that 


lim / Sk ■ Dkpdx = / Sod ■ Dua^ipdx. 

Jn\Ej Ja\Ej 

Consequently, since the graph is monotone, we observe, for arbitrary S e Rfym 
all nonnegative (p G L°^{n\Ej), that 

Os: lim/’ {~Sk-S^{;S)) :{bk-d)pdx. 

Jn\E^ 

Since p was arbitrary, we have that 

[Sao — S)) : {Duoo — (5) 5= 0 for all 6 G and a.e. x G n\Ej. 

According to Lemma this implies that 

(DMoo(a^), 5 ^ 00 ( 2 ;)) e A{x) for almost every x G ^\Ej. 

The assertion then follows from |ifj| ^ 0 as j ^ 00 . 

It remains to verify that Ofc 0 in measure as /c ^ 00 . We divide the proof into 
four steps. 

Step 1 : First, we introduce some preliminary facts concerning discrete Lipschitz 
truncations. For convenience we use the notation 


Ek := :iiUUk - «ao) = Uk- 4v“ao 6 V(0fc) 


and let be the sequence of Lipschitz-truncated finite element func¬ 

tions according to Proposition Recall from Lemma 1^ that Eu_^ 0 weakly 


in IFo’’'(0)“ 


i. e., we are exactly in the situation of Proposition 


16 


Although 


Ek e Yo{Gk), i-e., Ek is discretely divergence-free, this does not necessarily imply 
that Ekj G Yo{Qk) and thus we need to modify Ekj in order to be able to use it 
as a test function in (3.20). With the discrete Bogovskii operator := from 
Corollary [T^ we define 


■■= ^\diYEk,j) G YiGk). (6.8a) 

The ‘corrected’ function 

^k,j ■= Ekj — 'i’kj e Yo{Qk) (6.8b) 


is then discretely divergence-free. 

To this end we recall from Section 


We need to control the correction in a norm, 
that Q(C/fe) = spanlQ^,... ,Q%J for a 


3.2 


certain locally supported basis. Then, thanks to properties of the discrete Bogovskii 
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operator and Corollary [T^ we have that 


II'5' 


sup 
QeQ(t'fc) 


sup 


Q div Ek j da; Q div E^ j — div Ef^ da; 


IIQIIu 


= sup 

QeQ(^A;) 




Q-Hi^PiQi \suppQ'^cz{Ek,j=Ek} 

+ E 

sMppQ’^n,{Ek,ji=Ek]¥=0 


WQWr- 

f^p,Q^div(E>^’^ -Ek)dx 

IIQIIu 


PiQ\ d:vf{Ek,j - Ek) dx' 


sup 


E 


Q-Ei = i PiQ'i \s\ippQ^n{Ek^j^Ek}^0 


sup 


E 


Q=Y!^=iPiQ'l \snppQ^^n.{Ek,j¥=Ek}i=0 


< 


div-Efej-xn? 


{Efc,,■#£(,) 


sup 

Q=£,"\ PiQ’t 


PiQ'l; div{Ekp - Ek)dx\ 

K ) 

PiQ'^ div Ek,j dx\ 

) 

Ssupp Q^r\{Ek,j¥=Ek}9^0 P^Qi 


WQWr' 


^ C 


div£;fc,^Xn^^ 


^ C 


^Ek,j Xnif 




where Vo^ is the characteristic function of the set 


n 


{Et: j¥^E0 El \ E ^ Sk such that E cz {E^p ¥= -Efc}|. 


Note that in the penultimate step of the above estimate we have used norm equiv¬ 
alence on the reference space Pq from ( |3.1b I. In particular, we see by means of 
standard scaling arguments that for Q = Xii=i PiQi the norms 

/ / > \ '^ir 

(EI-°*IE|Q?IIu) and Q-IIQt, 


are equivalent with constants depending on the shape-regularity of Qn and Pq only. 
This directly implies the desired estimate. 


Observe that |r2|;| 


^ c 


c > 0; hence, 
that 


Qk 


E\ for all E e Qk, k eN, with a shape-dependent constant 
^ c\{Ekj 7^ Ek}\, and it follows from Proposition It |f) 


I3r\\^k,,\ 


l,r 


^ C 


^kjXn'f , 


^ c2-^> WEk 


(6.9) 


Moreover, we have from Proposition 16 and the continuity properties of *B" (see 
Corollary [T^ that 

weakly in IPq’®( 0)‘^ for all s e [1,cx;), (6.10a) 

^k,j,'^k,j 0 strongly in for all s e [l,oo), (6.10b) 


as k 


00 . 


Step 2: We shall prove (recall the last line of (6.6) for the definition of Uk) that 


limsup / jofel dx ^ c2 

n—>cc J {E^,=Eh,j} 


with a constant c > 0 independent of j. To see this we first observe that |a/c| = 
Ufc + 2a^ with the usual notation a^(x) = max{—afc(x),0}, x e 12. Therefore, we 
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have that 


limsup / 

/c—»-oo J {E^ = E^’^} 


\ak\ da; ^ limsup 


Ofe da; 


fe-’-oo a{£;"=£".3} 

+ 2 lim sup / 
fc—>-00 


( 6 . 11 ) 


Oj, da;. 


By choosing (j^u '■= Xsupp(a~) ^ in ( 6 . 61 , we observe that the latter term is 

zero. In order to bound the first term, we recall ( 6.81 and observe that 

f akdx=f (^fc - 5 *(-,Dmoo)) : - Dwoo) da; 

J{Ek = E^j} J{Ek = E;,j} 


+ / Sk -D^kjdx + / S'fciD^'fcjda; 

an an 

- f S*{-,Duao) : BEkj dx 

Ja 

+ [ {S*{-,Duoo)-Sk) :BEk,jdx 

a {Ek=ftEk,j} 

= Ifcj' + Ilfcj + Illfe.j + IVfcj + 


Thanks to (5.21 and (3.221 we have that 

f \Sk{-,DUk)- 


' {Ek = Ek,j} 

^ ||Sfe(-,DC/fe)-5*(-,D-u»)t, 


D^div^oo -BUo 


da; 


DDf^;';Moo -Buo 


0 , 


as k ^ CO. In order to estimate II^j- we recall that ^k,j e Vo(0fc) is discretely 
divergence-free, and we can therefore use it as a test function in (3.201 to deduce 
that 

llkj = -B[Uk, Uk, + / / ■ ^k,j da; ^ 0 as k ^ co. 

an 

Indeed, the second term vanishes thanks to (6.10al. The first term vanishes thanks 
to (6.3l and the weak convergence (6.10a) of ^kj- The term 111^^ can be bounded 
by means of (6.9); in particular, 

lim sup llllfcjl s: lim sup ||5(-,Di7fe)||,,, t < c2-^>, 

fc—»-GO k—*-CC 


where we have used (3.22). Proposition 16 implies that 

lim IVfc , = 0. 

fe—»-00 

Finally, by ( |3.22 1 and Proposition 16 we have that 

lim sup IVfcjl ^ lim sup (||5'*(-, Dmoo)!,., + l|S'fc(-, D17fc)||,,,) \\BEk,jX{Ek 9 ^Ekj}\\ 

fc—»-GO fc—»-GO 


^ c 


2-aA, 


In view of (6.11), this completes Step 2. 

Step 3: We prove, for any ■& e (0,1), that 

lim / lofcl’^ da; = 0, 

an 

which then implies the assertion —> 0 in measure as k —>■ co. 


( 6 . 12 ) 
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Using Holder’s inequality, we easily obtain that 

[ \akf dx = [ 

Jn Ji 


Ofcl’’ dx 


^ |U| 


l-'d I f 


lafcl’’ dx + 

J{- 

lofcl da; j + \ak\ da;^ \{Eu Ek^j}\ 


' {Ek¥=Ek,j} 


1 - 1 ? 


' {Ek=Ek,j} 


Thanks to (3.221, we have that (f^ |afc| da:)*^ is bounded uniformly in k and by 
Proposition 


16 we have that 


\\EX c 


A 


k,j 


221 r > 


where we have used that {Ek}keN is bounded in mQ’’'(f2)“ according to (3.221 and 
Assumption Consequently, from Step 2 we deduce that 

lim sup f dx ^ c 2“-^’^/’' -f 

A;—>00 Jq 


22^V(1-i9) ■ 


The left-hand side is independent of j and we can thus pass to the limit j —> oo. 
This proves ( |6.12[ ). □ 

6.3. Proof of Lemma |27[ Since rik ^ N as k ^ co, we may, w.l.o.g., assume that 
Uk = N for all fc e N. 

Step 1: We shall first prove that, in this case, we have that the (sub)sequences 
in Lemma [2^ do actually converge strongly, i.e.. 


"^00 

in Wj’*(L!)^ 


Pec 

in Lo(U), 

(6.13) 

= S^{Du^) 

in 



SkiDUk) = S^iDUk] 

To this end, we investigate 

Ofc := (S'^(DC/fe) - S^(Du^)) : D(i7fc - u^) > 0 

(compare with Assumption distinguishing two cases: r ^ and r > 
J, then we can conclude, as in the proof of Lemma 


If ^ ^ 

^ d+2^ 


that 


26 


in Section 


6.2 


Oii Jjak\^dx = {DU k)-S^iBu^)) : D(C/fe - Moo))’’ dx ^ 0, (6.14) 


where we have used that = 0 almost everywhere in U. Thus, recalling that S 
is strictly monotone, we obtain that 


N 


DUk Duo 


and 


^^(DJTfc) ^ S"(Dmoo) a.e. in U, 


tiV 


(6.15) 


at least for a subsequence of fc ^ oo. Since 1 ^ t < r and 1 sg t < r' (compare 
with (4.3al), we obtain with Lemma 24 and Vitali’s theorem that 


Uk 


inWo’'(U), and ^^(DC/fe)-> S'""(Dmoo) = S'oo in L*(U) 


N / 


\dxd 


Indeed, using Holder’s inequality, we obtain with (3.22) that 

||L)[/fc — DmqoIIj^^ ^ |w| rt ||DC/fc — Dttooll^^^ 

for all measurable w c U; i.e., {\DUk — E)itoo|*}fe6N is uniformly integrable and 
the claim follows from Vitali’s theorem and (6.15). The convergence of the stress 
sequence follows analogously and the claim S^'‘ (Dmoo) = Sco is a consequence of 
the uniqueness of limits. 
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If r > then t = r and i = f = r'; compare with (4.3b I. We deduce from 


(6.16) 


Lemma 24 and ( 3.23[ ) that 

0 sS limsup [ {S^{T)Uk) - S^{Du^)) : D{Uk - Mqo) dx 
A;—>00 Jq 

= limsup / 5''^(DJ7/c) : DC/fe — Soo : Dmqo da; 

A;— >00 Jq 

= limsup / f ■ Uk - Sao ■■ Dmoo dx= f ■ Uao - Sos ■■ Dmoo da; = 0. 

Ai^oo Jq Jq 

As before, thanks to the strict monotonicity of , we have that 

Uk Moo and S^{DUk) -> ^^(Dmoo) a.e. in n 

at least for a subsequence of fc ^ oo. Moreover, 

lim [ S^inUk) ■■'DUkdx = [ 5^(Dmoo) : Dmoo da;, 

Jn 

and thanks to (3.221, we have that S^(I)Uk) '■ DLffc is bounded in L^{Q), hence 
S^iDUk) : DLffc ^^(Dmoo) : Dttoo; compare with Lemma Recalling As¬ 
sumption we have 0 ^ m -f S^(DUk) ■ DC/^ almost everywhere in fl. Com¬ 
bining these properties, it follows from Lemma 0 that m + S^iBUk) ; DLTfc — 
TO -f S^{Duao) ■ Dttoo in L^{fl) and thus S^{DUk) ■ DC/fc ^ 5 '^(Dmoo) : Dttoo 
in L^(fl). Consequently, by the Dunford-Pettis theorem, : DUk}km is 

uniformly integrable. Thanks to the coercivity of , we have that {\Y)Uk rifeeN 
and {|S'^(DLrfc)| ’’ }fceN are uniformly integrable and hence we deduce from ( |6.16[ ), 
with Vitali’s theorem, that 

Uk Moo in and S^{DUk) S^o in 

It remains to prove the strong convergence of the pressure (sub)sequence {Pk}keN 
in L^(n). Thanks to (2.31, for A; e N there exists a, Vk s iS^Y with ||ufe||- = 1, 
such that 


otf boo - PkWi ^ / (Poo - Pk) div Vk dx. 
Ja 


Since {v/cjfceN is bounded in {^Y^ there exists a not relabelled weakly con¬ 
verging subsequence with weak limit v e ITg^’* (P^Y- Therefore, we deduce using the 
properties of (see Assumption ^ that 


ctr boo - Pk\\i^ / (Poo - Pk) div Vk dx 

Jn 

= / (Poo - Pk) div 3^^^Vk dx 4- / Poo div(t)fc - 3^^^Vk) dx. 

J Q J Q 

The second term vanishes thanks to Proposition For the first term, we have, 
thanks to Lemma and ( |3.23 1, that 

[ (Poo - Pk) div dx 

Jn 

= [ {Sao- S^iDUk)) : 3%Vk + (Bboo,Moo] - B[Uk.Uk\) ■ 3i^Vkdx 
Jn 

^0 as fc ^ 00 . 

Here we have used in the last step the strong convergence of {S^(I)Uk)}km and 
{Uk}k(Bf^ in Up)'^^'^ respectively WY*p)‘^, as well as that the latter result im¬ 
plies that B[Uk,Uk] B[ttoo,Moo] strongly in L* (H)"^. This completes the proof 

of ( |6T3t - 
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Step 2 : As a consequence of ( |6.13| we shall prove that 
n{Uk,Pk,S^{DUk)) strongly in (6.17) 

To this end we observe, for v e ITq’* that 


np‘^^{Uk,Pk,s^{BUk))-np^%u^,p^,Soo), v'^ 

= [ {S^{BUk)-Sao):Bv+{B[Uk,Uk]-B[Ua„Uao])-vdx 

Jq 

+ {Pk- Poo) divnda; 

Jn 

^ { ||5'^(Dt/fe) - Sao\\^+\\B[Uk,Uk] - B[Uao,U^]\\i+\\Pk-poo\\i}\\v\\,^i, . 


Hence, thanks to (2.1 1, (6.131, and the fact, that B[Uk, Uk] B[moo, Moo] strongly 
in L* (n) (see Step 1), this proves the assertion for The assertion for TZ"^ is 

an immediate consequence of (6.13). 


Step 3 : We use the techniques of |Sie09| to prove that 
Tl{Uao,Poo,S^{I)Uao)) = 0. 


(6.18) 


To this end, we first need to recall some results from |Sie 09]. 

For each x B the mesh-size hg^. {x) is monotonically decreasing and bounded 
from below by zero; hence, there exists an hoo e such that 


lim hg^ = hoo in (6.19) 

A;—>00 

compare e.g. with |Sie09l Lemma 3.2]. We next split the domain fl according to 
Gk ■■= f^Gi = {E e Gk : E e Gi for all i ^ k} and Gk ■= Gk\Gk, 

i^k 

i.e., setting := B.{Gk) and := 11(0)!), we have Gl = 11^ u 11°. It is proved in 
|Sie09l Corollarv 3.31 Icomoare also ill) that, in the limit, the mesh-size function 
hg vanishes on H^, i.e.. 


lim 


hgkXQO 


= 0 = lim 
A;—>00 


\hgkXni 


( 6 . 20 ) 


Here H), := and we have used the local quasi-uniformity of meshes for the 

latter limit. Since G^ c Gk Gk for any k ^ i, we have 


11° = 11(0°) = ll(0fc\e+)- 


Now, fix n e W^’‘'(H)^ n (H)'^ and q e 1T^’*'(H), with ||n ||2 p = 1 = ||g||^ 
We shall prove that 

(Tl{Uk.Pk,S^{BUk)). {v,q)) = (7^P°^(^7fc,Pfe,5'^(D^7fc)), 

+ <7^"(^7fc), (7-0^g> 


vanishes as A: ^ oo. Here we use the abbreviations := and := 0 q'‘. 
Then, ( 6.18| ) follows from (6.17) and the density of lT^’*(ll)‘^nVFQ’*(H)‘^ in 1 Tq’*( 11)‘^ 
and of (H) in L* (H). We shall estimate the two terms on the right-hand side 


























30 


C. KREUZER AND E. SULI 


separately. For the first one, we have with Corollary [20] that 


< Y, £^t^{Uk,Pk,S^{BUk);Ey 


EeQk 




s; E £t{Uk,Pk,S^{BUky,E) 
E^gk\Gt 

+ E £tk{Uk,Pk,s^{BUkyEY'' 
E^Gt 

<£iy{Uk,Pk.S^{BUuyQk\Qty'' 

^pde/yr fT\TT + 




Vt) '^'^div'^\\i>^u^k{E) 


\/v - vadi^t!| 




Wv - 


+ £^gf{Uk,Pk,s^(BUkygn i|v.- v.div-iit,,;,s(e+), 

where we have used Holder’s inequality and the finite overlapping of the U^^{E), 
E e Qk- In view of Lemma [M| and Corollary we obtain that 

£Pgy{Uk.Pk,S^{BUk);gkW) ^£^t^{Uk,Pk,S^{BUk)) < 1 . 

Recalling (|3.6|), we thus obtain from the monotonicity of the mesh-size function that 


'7^P^^(^7fc,Pfe, S^iBUk)), V - 

< \\hg^xntt + £tiUk,Pk,S^{BUkygtf\ 

A similar argument shows that 

<7^-(^7fe), g-a^q> < \\hgAa^,nt + £l{Uk;grf'' ■ 

Thanks to ( |6.20[ ), for e > 0 there exists an z e N such that 

(n{Uk,Pk,S^{BUk)), {v, q))<e + £g, {Uk, Pk, S^{BUk);g+), 
and it therefore remains to prove that 

£g,{Uk,Pk,S^{BUk);g+) ^0 asfc^o) (6.21) 

in order to deduce ( 6.18[ ). To this end, let 

£k ■.= £g,{Uk,Pk,S^{BUky,Ek) :=max{£g,{Uk,Pk,S^{BUkyE): EeMk). 

Then, by the stability estimate. Corollary and (|4.4b I we have that 


r,Ue(Ek) 


£k ^ kli^t-M^iEk) E l|f^fc|ll,t;WS(Efe) + \Pk\i-US{Ek) + 'if'^i,US{Ek) + 

+ ||div U k\f,us {Ek) 

S: ||Lffc — 'l‘co|li_t;iv6(Bfc) + k — '*^cx)||i_t;ZYS(Ej,) + \Pk — 'Pa)\i-;US{Ek) 

+ I|divi7fc||‘.j^g(£.^) 

+ (Ek) + ll“®lll.t;W0(Efc) + \Pa^li-U^ {Ek) + 'if\i,U^{Ek) + •• ^ifg(^Ek) 

The first line of this bound vanishes thanks to (6.131. Since Ek e £1^, we have that 
\Eky^‘^ ^ \hg^ lloo-no the remaining terms therefore vanish thanks to (6.201 and 
the observation that Ek e £1*1. Therefore, we conclude with that 

£g,{Uk,Pk,S^iBUk)-,gy) ^ #g+max[£g,{Uk,Pk,S^{BUkyE): Eeg+] 

^#gy9{£k)-^0 asfc^oo. 
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where we have used the continuity of g at zero and that c 
Combining these observations proves (6.181. 

Step 4: In this step, we shall prove that 

£g,{Uk,Pk,S^(DUk))^0 asfc^oo. 


Gk\Mk- 


To this end, we observe from (6.211 that it suffices to prove that 

Eg, {Uk,Pk, S^iDUk);gk\G+) - 0 as /c ^ cx) 

for some fixed i ^ 0. In view of Corollary 22 (6.171 and (6.181, it thus suffices to 
show that 

osc(C/fe,5'^(DC/fe);0A0+)-O 


as k 


00 . 


This is a consequence of the properties of the oscillation, (6.201, (6.13), and As¬ 
sumption noting that 

S^{BUk)-Tlg,S^{BUk) ^ ng,5ao-ng,5'^(Di7fc) 


i-,n° 


i-n° 


II He, 






hn? 


Observing that this readily implies that the estimator vanishes on the whole se¬ 
quence completes the proof. □ 


7. Graph approximation 


In this section we shall discuss the approximation of certain typical maximal 
monotone graphs satisfying Assumption Admittedly, for particular problems 
the approximations suggested here might not always represent the best possible 
choices, and in the context of discrete nonlinear solvers, such as Newton’s method, 
properties of the smoothness of the approximation may become important as well. 
We believe however that the following examples provide a reasonable guideline for 
constructing graph approximations with properties that are required in applications. 


7.1. Discontinuous stresses. Typical examples of discontinuous dependence of 
the stress on the shear rate are Bingham or Herschel-Bulckeley fluids. In this case, 
the fluid behaves like a rigid body when the shear stress is below a certain critical 
value S* and like a Navier-Stokes, respectively power-law, fluid otherwise; compare 
with Figure [2 To be more precise, for some yield stress cr ^ 0, we have 
l-S'l ^ cr D = 0, 

|5|>a ^ 5 = a^ + 2K|£»P)£>; 

where u > 0 denotes the viscosity u > 0; see |DL76| . A selection of the correspond¬ 
ing maximal monotone graph is given, for example, by 


S*{D) :=S^\D\) 


D 

JD\^ 


with S* (D) := 


|o. 

\a + 2v{D)D, 


if £» = 0 
otherwise. 


(7.2a) 


For the sake of simplicity of presentation, we restrict ourselves in the following to 
u > 0 being a constant. However, we emphasize, that the approximation techniques 
presented below can be generalized to more complex relations such as, for example. 


S*{D) 




if D <5 
otherwise. 


S:{D)=c,{K^ + D^y-^D, (7.2b) 


for D ^ 0. Here 6^0 and ci, C 2 , ki, K 2 5= 0, qi > 1, 52 = t, such that S'i((5) ^ S 2 {S). 

We denote the maximal monotone graph containing {{D, S*{D )): < > 0} by a and 
observe that {S,(t) e A if and only if (|^|, |(t|) e a. Therefore, the approximation of 
the monotone graph reduces to approximating the univariate function S* by some 
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smooth 5” : Kj —> M(|. The explicit smooth approximation of S* is the obtained 
setting 

S-{D) := S-{\D\)^^ for all D e (7.3) 




Figure 1. Bingham fluid (left) and schematic approximation of a 
more complex law (right). 


A simple approximation. A simple approach to approximating S* in (7.2a) is to 
use the following smooth explicit law (cf. [HMSTl^ l: 

S'^{D) := (2v A ——)T>, where Dt := 'J+ r^. 


First assume that a > S'^{D); then, {0,S'^{D)) e a according to (7.1). If I? ^ r, 
then \S'^{D) — S'^{D)\'^ + \D — Op = \D\'^ ^ t^. Otherwise, we have 


{2i^+—)D 


and hence 


(TT^ > 2vD Dr{D + Dr) > , 




\S%D) - S^{D)\^ +\D- Op = |Z7p ^ r 


Assume now that 0 ^ u < S'^{D); then, D ^ t implies 


2j/£» +cr^ = S'^(L>) > cr —D> 

Dr <J 


^ 1 

Dr{D + Dr) ^ 4' 


In other words this case can occur only for ‘large ’t ^ If Z? > r then we have 
similarly 


< 2i'D Dr{D + Dr) ^ 8vD^ 


\8iyJ 


N 1/3 


Therefore, we obtain 

\S^{D)-S^D)\ = a 


Dr-D 

Dr 


< 


Dr{D + Dr) D‘^ 


^ < 41.2/3^1/3^2/3^ 


Combining the above cases shows the validity of Assumption 28 with r = -, for 
example. The verification of Assumption is left to the reader. 

Approximation by mollification. We can extend S* to an odd function on the 
whole real axis by setting 

S*{D) ■=-S*{-D) forZ)<0. 

Then, for n e N, we define an approximation of S* by 


/ oo 

5'*(s) 77"(s-t) ds, 

-00 
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with = nr](nt); here rj e is a nonnegative even function with support 

(—1,1) such that ds = 1. Consequently, the function S'” e C(K) is odd and 

thus S”(0) = 0. 

For nonnegative D e Kj we have by the monotonicity of S* and the definition 
of S”, that there exists some 0 ^ D* e {D — + d), such that {D*, S'^{D)) e a. 

Therefore, we have 

l^np) _ 5'n(£))|F + |£) _ ^ 0 + ^, 


and consequently 


S^(5,S”(^))^-^0 asn 


00 


for all S 6 


< d 
^"sym- 


This shows that Assumption 28 holds. Moreover, tpn satisfies 


Assumption]^ compare e.g. with |BGMSG09l IGMS^ IGZG07] . 


7.2. Monotone graph with plateaus. Similarly to (7.21, we consider a maximal 

i|^, but now assume that S* : 

is continuous with 


monotone graph with selection S*{D) = S'*(|Z)|)| 


^0 


^0 


S%D) 


< (T = const., 


if D < Si 
if 6 i D < 62 

else. 


with 

S*iD) =Ci{K‘f + D^)‘^D, i=l,2. 

Here ci, C 2 , ki, ^2 > 0, qi > 1, <72 = r, such that Si (5*) = S'J = compare 

with Figure (left). In this case, we are basically in the same situation as in 




Figure 2. Graph with plateau (left) and plateau and jump (right). 


Section [7.1| with interchanged roles of S and D. Therefore, using the approximation 
techniques of Section |7.1| we can construct an approximation of the monotone 
graph where the shear rate depends explicitly on the shear stress. However, in a 
practical numerical method this relation typically has to be inverted, which may 
cause additional computational difficulties. 

Another approach is to use an approximation of the form 


S^{D) := 


< 82 ( 0 ) 
S^iD) 


if S^{D) < a-^ 
if SiiD)>a+^ 
otherwise. 


where is the linear interpolant between cr — d and cr + d with corresponding 
values for D. 
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Combined with an approximation strategy as in Section [7.1| this procedure can 
also applied to cases where jumps and plateaus are both present; compare with 
Figure]^ (right). 

Remark 32. The arguments of §^7.1\ and \ 1 . ^ can obviously be extended to finitely 
many jumps/plateaus and in particular cases even to countably many jumps/plateaus. 
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Appendix A. Stability of the a posteriori Estimator 
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Referring to Remark 
belong to Uoc(^) when 1 


we shall illustrate that divS"(DC/g) does not need to 
< r < 2. In fact, this problem already appears in the 


a posteriori analysis of quadratic finite element approximations of the r-Laplacian. 


To see this, we consider the following problem: for / e L’' (fl), find u e WQ’^(r2), 
such that 


■ div (iVur-^Vw) = / e 


For Q e G, £ ^ 2, let Ug e V(^) = {V e C(fl) : V\e e for all E e Q}. Then, for 
the residual 


7^(^7) = / + div(|VC/r-2Vf7) e (G) 
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and r ^ r', we have the bound 


divIVC/ar-^VC/a + < C ^ \\hg{-div\yUgr^yUg + 

Eeg 


E ||[[/^f'|Vt/gr-"VC/g]] 


Eeg 


r 

r,8E 


:= E S{U,g,Ey =:£{U,gf; 

Eeg 

compare e.g. with |BS08] . In fact, it can be shown as in |Verl3l Chapter 3.8], that 

\\niU)\\_^ .<£iU,g) < ||7^(C7)||_^__ + osc(t/,0). 

However, it is not clear, whether the element residual 

E ||/ie(-div|VC/gr"Vt/g + /)||;^ 

Eeg 

(and thus £{U, g) and osc{U, g)) is bounded. In fact, elementary calculations show 
that on an element Eeg, we have 


divdVC/r-^VC/) 

= E 


vt/ 

|W|2 


(VC/).) + (r - 2)\vur^j^{vuy] 


d^c/ 

dxf ' 


8 ^U 

{dxi)'^ 


co,E 


As a consequence of U\e e P^, we have that 

T j 

and may have different zeros, we need to 

f dx < ao. 

Je 

Since Vt/ e P^_i, for this to hold we need to have that 

(C-l)(r-2)f >-1 ^ £<l + 


|VC/|2' 
is bounded. Since |VC/| 


1 


(A.l) 


r(2 - r) 

where we have used that 1 < r < 2. This puts a restriction on the polynomial 
degree of the finite element space. 

We consider some special cases: 

• d = 1 and C=2, l<r<2 and r = r'. For simplicity, we assume that 
E = [0,1] and that U = x'^. Thus VC/ = U' = 2x, div VC/ = U" = 2, and 
we have that 

f |div(|VC/r-2vc/)r'da: = f I div(|2a;r-22a;)r'dx = 


^2C-l)r' r|div(x’-l)|"'dx 
Jo 

= da;. 

Jo 


This term is finite if and only if 

(r — 2)r' > —1 
r(r - 2) > 1 - r 
r^ — r—1>0 
, A - V5 1 +V5. 

^ —f;—,—f;—]• 


2 


2 
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This restriction applies to any dimension, since one could always choose 
U = [xi, 0,..., 0]* as a quadratic function on the standard d-simplex. 

1 < r < d = 2 and f = This case is relevant for power-law type 

models when r is close to 1, or, more precisely, when r ^ respectively 


^ compare with (2.1l. Then, (A.l) becomes 




Since r > 1, in this case the residual is in general not in even for 

quadratic elements {i = 2). 
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